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We study the interior structure of the Einstein-Yang-Mills-Dilaton black holes as a func- 
tion of the dilaton coupling constant 7 e [0,1]. For 7 7^ the solutions have no internal 
Cauchy horizons and the field amplitudes follow a power law behavior near the singular- 
ity. As 7 decreases, the solutions develop more and more oscillation cycles in the interior 

' region, whose number becomes infinite in the limit 7^0. 
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I. The much discussed mass inflation scenario |lj provides a picture of what happens to 
^ , the unstable Cauchy horizons of the Kerr-Newman family when the back reaction of the 

blue-shifted perturbations is taken into account. In this scenario the blue-shifted influx, 
due to the inavoidable gravitational wave tail of a realistic gravitational collapse, causes 
the local mass function to blow up exponentially near the Cauchy horizon. Thereby, a 
singularity is produced along the Cauchy horizon, which leads to infinite tidal forces on 
a free-falling object. This singularity is, however, lightlike and mild, accompanied with 
a divergence of the conformal curvature, but with bounded shear and expansion. The 
genericity of the mass inflation picture has been established in a series of analytical and 
^ ■ numerical studies of a class of vacuum and electrovacuum solutions [[|. 

What really happens very close to the singularity along the horizon cannot be de- 
cided by classical physics. This is presumably not even a physically relevant question, 
because the black hole will have evaporated or merged with other black holes in a cos- 
(3J[)! mological big crunch. It may, nevertheless, be of some interest to study the interior 

black hole solutions for classical field-theoretical matter models which are important in 
high-energy particle physics. 

Such investigations were initiated in ref. [El, where it was shown that for pure 
Yang-Mills fields a new type of infinitely oscillating behavior with exponentially growing 
amplitude is developed. It is natural to ask, whether this cyclic repetition with associated 
mass inflation is generic when other matter fields are included. In J3|, [||] it was shown 
that the interior solution changes qualitatively when a Higgs field is included. 

In the present contribution we discuss the interior structure of black holes for the 
Einstein-Yang-Mills-Dilaton (EYMD) system, considering the dilaton coupling constant, 
7, as a bifurcation parameter of the associated dynamical system. It turns out that for 
7 = 1 the behavior is drastically different from that of the EYM system (7 = 0), in 
that there are no oscillations; this phenomenon has also been observed in f|. All fields 
evolve quite monotonically when the radial Schwarzschild coordinate r (internal time) 
approaches 0. This has prompted us to study the change of the phase portrait as 7 
decreases from 1 to 0. Our analytical and numerical results show quite convincingly, 
that with decreasing 7 more and more cycles are developing, and that the limit 7^0 
is quite singular. For 7 7^ the interior solutions have no Cauchy horizons. We believe 
that this is true quite generically for a large class of nonlinear field-theoretical matter 
models. Further work on this issue is in progress. 
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II. The black hole solutions we shall be considering arise within the context of the SU(2) 
EYMD theory, whose action reads in standard notations 

S = J (~R+^d^d^ + ^e 2 ^trF, v F^ ^d A x. (1) 

The dilaton coupling 7 is assumed to have values in the interval [0, 1]. We are especially 
interested in the limiting behavior for 7 — - > 0. When 7 vanishes, the dilaton decouples, 
in which case one can put $ = 0. In the static spherically symmetric case the metric is 
given by 

ds 2 = S 2 Ndt 2 - ^ - r 2 (d$ 2 + sin 2 dip 2 ) (2) 
and the gauge potential A for a purely magnetic YM field can be parametrized as 

A = w(r)(-T 2 d9 + Ti sin 9 dp) + T 3 cos 9 dip, (3) 

where the group generators T a are chosen as r a /2i (r a =Pauli matrices). The functions 
w, S, N and the dilaton $ depend only on r. It will be sometimes convenient to express 
N in terms of the mass function m as N = 1 — 2m/ r. The field equations, corresponding 
to (|ip read 

(rN)' + r 2 N& 2 + U = 1, 
(SNr 2 &)' = jSU, 

r 2 (NSe 2 ^w')' = Se 2 ^w(w 2 -l) 1 

S' = S(r(P' 2 + 2e 2 ^w' 2 /r), (4) 

where U = 2e 27 * (Nw' 2 + (w 2 — l) 2 /2r 2 ). These equations are invariant under the scale 
transformations 

r i-> e A r, N ^ N, S h-> e~ x S, $ h-> $ + -, w i-> w. (5) 

7 

The function S can be eliminated from the system When N, $ and w are known, 
S can be expressed as 

S = exp (- J r °°(r& 2 + ^e 2r V 2 )c^) . (6) 

For any 7 e [0, 1], equations @ are known to possess a family of black hole solutions 
0, [@. Close to the event horizon they behave as follows: 

N=(l-V h )x + 0(x 2 ), $ = $, + -^-x + 0(x 2 ), w = Wh + ^l^l x + 0(x 2 ), 

U — Vh) U — Vh) 

(7) 

where x = (r — r^/r^, = e 27 * h (u^ — l) 2 /?" 2 ., and the scaling symmetry (H) can be 
used to set $h = 0. Note that at the event horizon one has N'(rh) > 0. The asymptotic 
behavior at infinity is described by 

ir „ 2M D 2 _ ^ D 1 / c\ 1 

iV=l + — + 0(-), $ = $ 00 -- + 0(-), w = ±(l-~)+0(- 
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Figure 1: The interior solutions for the n — 1, — 2 EYMD black holes. On the left: 
for 7 = the metric functions m and S can be approximated by step-functions whose 
amplitudes grow exponentially as r — ► 0. For 7 = 1, m and S are described asymptotically 
by linear functions - on a logarithmic scale. The picture on the right shows how the steps 
transform into straight lines as 7 increases. 



where M is the ADM mass and D is the dilaton charge. The black hole solutions 
with these asymptotics are numerically known in the whole interval ^ < r < 00 0, 
||. Their behavior is qualitatively similar for all values of 7: the functions N and $ 
are monotone in the exterior region, interpolating between the boundary values given 
by Eqs. (^), (d), whereas w oscillates within finite bounds. The solutions form a 2- 
parameter family labeled by r^, and n, the number of nodes of w in [r^, 00). In what 
follows we shall consider the extension of these solutions into the interior region of the 
black hole, < r < r^. It will turn out that the behavior of the solutions in the interior 
region changes dramatically with varying 7. 

III. When extending the black hole solutions into the interior region, the following 
phenomenon is observed: For 7 = (the dilaton is absent) the interior solution is 
characterized by violent oscillations of the metric coefficients and the gauge field strength 
0. The amplitude of these oscillations tends to infinity as the system approaches the 
singularity. For 7 = 1 the behavior is rather different, without any peculiar variations 
(see Fig.l). It is interesting to study more closely what happens when 7 decreases form 
1 to 0. Our analytical and numerical investigations lead to the following qualitative 
picture: For small values of 7 the solution in the vicinity of the horizon follows closely 
that for 7 = 0. The smaller the 7, the more oscillation cycles the solution exibits. 
However, as long as 7 7^ 0, the functions assume a power law behavior in the vicinity of 
r = 0: 

N oc -N ir ~ {1+p2) , $ oc $1 + p ln(r ) , w (x w - &r 2(1 ~ 7p) , (9) 
with b,Ni,p > 0, and the values of p are restricted by 



7 2 + l-7<p<-. (10) 

7 

In order to gain some qualitative understanding of this phenomenon, it is helpful 
to introduce the new variables 

x = w'e^, y = - (w2 ~ 1)2 e 2yI> , z = r$'. (11) 

r 2 N 
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Figure 2: On the left,: solutions of the dynamical system (13) forj = and 7 = 0.00005. 
For 7 = the trajectory is bound to stay in the z = plane always spiralling around 
the center, whereas for 7 = 0.00005 it leaves the plane and converges to the z-axis. On 
the right: the solutions for 7 = 0.005 and 7 = 0.020. For 7 = 0.005, one additional 
oscillation cycle is performed. 



In addition, the following approximations in the differential equations are numerically 
justified for the interior solutions: 

• w ~ const, (but w' is kept), 

• e 2 ^(w 2 - l) 2 > r 2 , 

• the term w(w 2 — l)/r 2 in Eqs.(H) can be neglected. 

This yields, as an approximate first integral of the field equations, 

SNe 2 ^w' ~ const. (12) 

Eqs. (H) can then be truncated to the following dynamical system: 

x = x(y + 72 — 1), 

y = -y(2x 2 -y-l + z 2 + 2 1 z) 1 

z = - 7 (2a; 2 - y) + yz, (13) 

where a dot stands for differentiation with respect to t = — ln(r). We are interested in 
the behavior of the solutions as t — » 00. First of all, we note that the planes {x = 0} 
and {y = 0} are invariant sets and hence divide the phase space into four regions. Since 
N is negative and the equations are symmetric under x 1— > —x, we restrict ourselves to 
the invariant region x > 0, y > 0. The critical points in this region are: 

1. x = a/3, y = aft, z = —07, where a = 1/(1 + 27 2 ) and (3 = 1 + 7 2 . 

2. x — 0, y — 0, z — p, where p is arbitrary. 
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The eigenvalues corresponding to the first critical point are Ai = a/3, A 2i3 = a(3(l ± 
+ 327)/2, hence it is an attractive center in the limit t — » — oo. The solutions 
starting at this point spiral outwards when t increases. For 7 = one has z = and 
Eqs.fljjD reduce to the 2- dimensional system 



x = x(y- 1), 

2 / = l/(-2x 2 + i/+l). (14) 

Now, for this system one can show that the spiraling motion around the focal point can 
never stop. The trajectories are bound to spiral around the center (1, 1) for all t. At 
the same time, they have to remain in the region x > and y > 0. As a result, after 
each revolution the trajectories come closer and closer to the second critical point, (0, 0), 
which corresponds to the spacetime singularity, but can never be reached. This explains 
qualitatively the unbounded oscillations for the pure EYM case ||. It is interesting 
to observe that the critical point around which the spiralling occures has no physical 
meaning by itself. Indeed, it can only be reached for t — > —00, which corresponds to 
r — ► +00, contradicting the assumption r < r^. 

On the other hand, for 7 7^ the solutions can no longer stay in the plane {z = 0} 
for all t. It turns out that the appearance of the additional degree of freedom completely 
changes the behavior of the system. When 7 is small the trajectory can stay for a 
long time in the vicinity of the plane spiralling outwards. The smaller 7, the smaller 
z and more revolutions are executed by the trajectory. As long as p = z is small, the 
second critical point is repulsive, since the corresponding eigenvalues are Ai = 1 — jp, 
A2 — P 2 + — 1, and A3 = 0. For small p one has Ai > > A 2 . However, when 



p = z becomes large enough to fulfill the condition (10), one will have Ai > 0, A 2 > 0, 
and so the critical point (0,0, p) will become attractive. One can also show that it is 
stable, despite the zero eigenvalue. As a result, having performed only a finite number 
of oscillations, the trajectory will eventually be attracted and end up at the singularity. 
The linearized solution near the critical point reproduces the power law behavior given 
by Eq. (§. 

IV. The interior solutions decribed above do not have inner horizons. It is natural to 
ask whether this property persists for all values of n and r^. For 7 = some of the the 
solutions are known to display Cauchy horizons. Specifically, for each n > 1 there is one 
distinguished black hole solution with radius r^(n), which has an inner horizon at some 
r = r_(n) 0. It turns out that for 7 7^ this can never happen. The proof of this 
statement is fairly simple and goes as follows: When 7 7^ 0, the scaling symmetry fl5|) 
implies the existence of the conserved Noether current 

J = r 2 (^SN' + NS' - NS ^j ■ ( 15 ) 

The conservation condition J' = can be straightforwardly verified with the help of Eqs. 
(|4|). Taking the limit r — * r^, we obtain 

2 

J = jS h N\r h ) > 0, (16) 

since S(r) is everywhere positive, and N' > at the event horizon. Now, assuming that 
an inner horizon exists at some r_ < r^, such that N(r_) = 0, N'(r_) < 0, and taking 
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the limit r -> r_, we arrive at 

2 

J = y5(r_)iV(r_) < 0, (17) 

which contradicts ([16]). The addition of a dilaton field thus eliminates the Cauchy hori- 
zons. We believe that the same is true for a large class of non-linear matter models. 

A more detailed account of this work is given in the diploma thesis by OS 0. 
MSV thanks the organizers of "The Internal Structure of Black Holes and Spacetime 
Singularities" workshop at Haifa for an interesting meeting and the generous hospitality. 
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